Usando un solver con condiciones de borde

Búsqueda por fuerza bruta variando los parámetros $c$ y $d$

C:\Users\rafae\AppData\Local\Temp/ipykernel_18660/1258265208.py:2: RuntimeWarning:

invalid value encountered in sqrt